Center-of-Mass Properties of the Exciton in Quantum Wells 



A. Siarkos, E. Runge, and R. Zimmermann 
Institut fur Physik, Humboldt- Universitat zu Berlin, Hausvogteiplatz 5-7, 10117 Berlin, FRG 

(February 1, 2008) 

We present high-quality numerical calculations of the exciton center-of-mass dispersion for 
GaAs/AlGaAs quantum wells of widths in the range 2 — 20 nm. The k -p -coupling of the heavy- and 
light-hole bands is fully taken into account. An optimized center-of-mass transformation enhances 
numerical convergence. We derive an easy-to-use semi-analytical expression for the exciton ground- 
state mass from an ansatz for the exciton wavefunction at finite momentum. It is checked against 
the numerical results and found to give very good results. We also show multiband calculations 
of the exciton groundstate dispersion using a finite-differences scheme in real space, which can be 
applied to rather general heterostructures. 



I. INTRODUCTION 



Excitons dominate the optical properties of low-dimensional semiconductor heterostructures such as quantum wells 
(QW) and quantum wires. The relative motion of the constituent particles and their center-of-mass (COM) motion 
determine different characteristics of the optical spectra and exciton kinetics. 

The exciton relative motion in QW is well studied and understood. The confinement of the carriers along one 
or two spatial directions into regions comparable to or smaller than the bulk exciton size enhances the effect of the 
electron-hole Coulomb interaction. This results in larger binding energies and oscillator strengths and in an increased 
stability compared to bulk excitons. Therefore, excitons are observed even at room temperature in these structures. 
The effect of the reduced dimensionality is as a rule much larger on the exciton groundstate than on the excited states. 

Details of the excitonic optical spectra of QW related to the COM motion like, e.g., inhomogeneous broadening and 
Stokes shift between photoluminescence (PL) and absorption are frequently used for structure characterization. These 
features are influenced by exciton localization and diffusion in the presence of interface or alloy disorder [jlj. Optical 
spectra and their temporal evolution are determined by the exciton formation processes [Q] and the subsequent energy 
and spin Q| relaxation dynamics. Spatially resolved spectroscopy techniques like micro PL and near- field scanning 
optical microscopy allow direct observation of exciton COM quantization in local potential minima |Q. All these 
phenomena are intimately related to the exciton COM properties whereby different energy and COM momentum 
regions of the exciton dispersion are probed in different processes. 

In many heterostructure systems of interest like, e.g., GaAs/AlAs, InGaAs/InP and ZnCdSe/ZnSeS, the exciton 
can be described in the effective-mass approximation (Wannier exciton) due to its small bulk binding energy (e.g., 
4 meV for GaAs). In this approximation, the degeneracy of the valence bands at the center of the Brillouin zone 
for materials of cubic or zinc-blende symmetry was first taken into account by Dresselhaus He also pointed out 
the absence of a well-defined COM transformation due to this degeneracy. Altarelli and Lipari j(| calculated the 
exciton COM dispersion for direct- and indirect-gap bulk semiconductors. They demonstrated that the ambiguity in 
the choice of the COM transformation can be used to achieve formal simplicity or optimal numerical convergence. 
For bulk GaAs, where the heavy- to light-hole mass ratio is large, the exciton dispersions are found to be strongly 
anharmonic and show avoided crossings between different branches. 

In semiconductor QW, the broken translational symmetry in the growth direction leads to the splitting of heavy 
and light holes at the T-point, and subsequent formation of heavy- and light-hole excitons. Due to the large hole- 
to-electron mass ratio, the influence of the valence-band mixing on the COM motion is greater than on the relative 
motion. The exciton dispersions are, thus, strongly non-parabolic. Direct consequences of the exciton dispersion 
anharmonicity in QW like slow indirect excitonic transitions due to camel-back shaped dispersions have been 
experimentally observed ||. 

The multiband exciton (i.e., with the full coupling of heavy- and light-hole bands taken into account) has been 
theoretically and numerically thoroughly investigated at vanishing COM momentum Q [p|-|l6||. The numerical effort 
for such calculations remains reasonable due to the high symmetry of this point. In contrast, only a few publications on 
multiband calculations of exciton COM dispersions in QW exist ( [17-19 ) since these are very demanding. Methods 
for improving the numerical accuracy and reducing the effort of such calculations are clearly necessary. Particularly 
useful would be an easy-to-use approach that gives the main features of the exciton dispersion with at least moderate 
accuracy. 
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The main focus of the present work lies on the exciton groundstate dispersion and its properties. A secondary goal 
is to study the feasibility of numerical methods for calculating the exciton groundstate dispersion in more complicated 
structures like V-groove quantum wires (poj , Excitons in GaAs/Al .3Gao.7As symmetric QW are considered (section 

f). Two different methods for the dispersion calculation are used: (i) the extension for Q ^ o f the well-known 
3[ expansion of the exciton in the product spa ce of electron and hole subband states (section I II Ah , and (ii) a 
finite-differences scheme in real space (section [IB ) with a groundstate-adapted Coulomb discretization |19| . Method 
(i) gives high-quality numerical results but is not feasible, e.g., for quantum wires of complicated geometry. With 
method (ii) the groundstate exciton dispersion in a V-groove quantum wire is tractable. Its convergence properties 
are checked here against the results of (i). Results for V-groove quantum wires will be presented elsewhere. 

Improving own previous results Jl9|] , we address in detail the problem of the choice of the COM transformation and 
introduce an optimized, groundstate adapted COM transformation that greatly enhances the numerical accuracy and 
stability of our results (section [II). 

Quasi as a by-product, a semi- analytical expression for the average exciton groundstate mass suitable for exciton 
localization problems is derived (section IV). This expression is of great practical importance since it gives for not too 
wide QW reliable mass values. The only necessary ingredients are the lowest subband dispersion and a good estimate 
for the groundstate Bohr radius. 

Finall y, we discuss the results of the exciton dispersion calculations in momentum- (section V A ) and real-space 
(section VB). The results of our semi-analytical expression for the average exciton groundsta te m ass are compared 
to various other mass expressions as well as with the numerical dispersions separately (section V C). 



II. THEORETICAL MODEL 



We consider the well-studied system of direct Wannier excitons in a single symmetric GaAs/Alo.3Gao.7As type I 
quantum well grown in (100) direction. Many aspects of the presented results can be effortlessly extended to Wannier 
excitons in other, more general heterostructures. 

In the envelope function approximation, the Wannier exciton is described by the Hamilton operator 

H = H e {r e ) + H h {f h ) + Vcoui(r e - r h ) . (1) 

h describe the material-dependent bandstructure of the respective particles in the vicinity of the T-point, and 
Vcoui stands for the attractive Coulomb interaction. We choose the coordinate system as usual with the z-axis 
in growth direction (100); r e — (x e ,y e , z e ) , rh — (xh,yh, Zh) denote the space coordinates of electron and hole, 
respectively. For the materials involved, the conduction band is to a good approximation parabolic; anharmonicities 
in the conduction band arise mainly through the interaction with the light and split-off valence bands which is small 
due to the relatively large band gaps. The valence band is adequately described by the Luttinger Hamiltonian pl| ] 
in the axial approximation ||,^2j , which takes into account explicitly the coupling of the heavy- and light-hole bands 
(rg) but suppresses warping. The coupling to the split-off (r?) band can be safely neglected for subband states 
with energies up to approximately 50 meV from the band edge because of the relatively large energy separation. We 
neglect the effect of the different dielectric constants (no image charge effects) [|l2| , ^5| , p3[ , and also all effects that lead 
to a small spin-splitting like lack of inversion symmetry of the bulk material |14[| or the interfaces |24| as well as the 
exchange part of the Coulomb interaction p7| . The electron spin is irrelevant and will be given a fixed value of +1/2 
in the present work. The quantization axis of the electron spin and of the hole angular momentum J is taken along 
the growth direction, and we use for the valence band edge states the same convention as in Ref. [^5| . 

The Hamilton operator H, Eq. (|l|), acts within these approximations on a four-component envelope function in 
the product basis of the conduction and valence band edge states {|| m,j) v |^ + ^) c }, where the hole-spin projection 
attains values of uij = +|, +i, — |, — |: 

H e = (-&Sdl + 91 + dl) +V c )l, (2) 
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I is the 4x4 unity matrix. The material parameters 7i(z/i), ^(zh), 73(^/1) a s well as the offsets V c (z e ) are 

piecewise constant functions of z e , z^. To ensure that the kinetic operators remain Hermitian in the presence of 
interfaces, we use the symmetric substitutions 

jdi ->• (<9i7 + 7<9j)/2, 7 9?- -> (9»79,- +dj'yd i )/2, i = x,y,z. (5) 

The in-plane COM momentum Q = —ih(\7 eil + V/j,,) is a constant of motion because the interaction term (^) 
depends only on the relative distance of the two particles. Reflection with respect to the central xy plane, a xy , is 
also a symmetry element for symmetric QW. Consequently, the exciton can be characterized by the parity P = ±1. 
Then, the wavefunction factorizes into 

3.) = £ **/-(z e ,z h)/ 5) I § mj) I i|) c , (6) 

mj 

where p = r e|| — 7%,, is the in-plane particle distance, R the COM space coordinate canonically conjugate to Q, and a 
stands for the remaining quantum numbers related to the relative motion of the exciton. 

The COM space coordinate R in is not unambiguously defined because of the anharmonic dispersions of the 
constituent particles [^J. The COM transformation must be linear in order to preserve the canonical commutation 
relations of space and momentum operators, and it has in general the form 

i? = /3r e|| +(1-/3) f h|| , fc = -ifc ((1 - /3) V e|| -pV h]] ) . (7) 

In the parabolic case, the free parameter (3 is taken as the scalar 

ftparab (8) 

m e + mh 

in order that relative and COM motion completely decouple. For bulk excitons, j3 has been considered in the literature 
as a scalar, a tensor in real space or even a spinor (26). We will return later on to the problem of an appropriate 
choice for the COM coordinate R. 

Taking into account the electron spin degeneracy, each exciton state is at least fourfold degenerate. It can be shown 
in a similar way as has been done for the hole subband states in Ref. |25| that the operator R-^T, with the rotation 
by 7r about the z axis, R^, and time-reversal, T, transforms between the degenerate states of different parity and 
opposite electron spin. If one combines this operator with the Pauli matrix a y , which flips only the electron spin, we 
have apart from an overall phase: 



vp* ~ Pa (f e , f h ) = (a% i? w r) p*(? e , r h ) ( 9 ) 

= P -^T £*£ ; / Wfc,-p)ll ~mj) v |i+|) c . 



Comparing m) with (O), we find: 



*%- Pa (z e ,z h ,ft) = ^J P f{z e ,z h ,-p) . (10) 



That is, the (degenerate) state of reversed parity is obtained by inverting the order of the spin components of the 
exciton envelope, complex conjugating, and changing the sign of the in-plane relative coordinate. Thus changing 
the multiband exciton parity with fixed electron spin in symmetric QW corresponds to flipping the hole spin in the 
single-band exciton case. In the axial approximation and for Q = 0, the different angular momentum components 
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decouple p2| , and changing the sign of p in Eq. (|To|) just changes the sign of two spin components leaving the other 
two unchanged; this holds no longer at Q ^ 0. 

We have solved the eigenvalue problem ([!]) in two ways which will be discussed in turn: (i) in fc-space, expanding 
Eq. (|(|) in the product space of the electron and hole subband states |h|, and (ii) in real space, using a finite- 
differences scheme. The first method gives very accurate results and is used to reveal the main features of the exciton 
dispersion. The second method is only suitable for the groundstate dispersion but promises to be feasible for more 
general structures. It is validated by comparing its results with the ones from the first method. 

A. Solution in fc-space 

As a first step, we calculate the single-particle subband states and their dispersions 

H e \n e k e \ zb 2) = £n e (^-e) |^e k e \ ^ 2) 

H h \n h k h ; p h ) = £ nh (kh) \n h k h ; p h ) (11) 
using a transfer-matrix method as in [ pT^ . The respective solutions in the axial approximation are of the form 

ik e -r 



n, l^^h = ^^e^ e ^ nc (ze)\^ ±i)c- (121 



ik h -r h 



n h k h ; Ph = ±l) = £ e im ^^f tphtkh (z h )\lmj) v 
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In Eqs. ( p"l| , |l2| ) n e ,ft denote the subband indices, fc e ,ft = (k e ,h, @e,h) the respective in-plane wavevectors in polar 
coordinates and p^ the hole parity under a xy [ B5J . 

In a second step, the exciton wavefunction for a given COM momentum Q is expanded into 

&2' Pa (f e , fh) =E / * ^InM K ^ \ + |>c \n h k h \ Ph)v , (13) 
n e n h 

with subband states of the two particles combined in such a way that the resulting exciton state has the required 
parity P and total momentum Q: 

k e = k + P$, 4 = fe-(l-/3)<9, P=p h -(-l) n ° +1 . (14) 

The last relation reflects that the conduction subband envelopes are even (odd) for odd (even) subband index. Fixing 
exciton parity P and electron spin eliminates any degeneracy at Q ^ 0. 

With the expansion ( |l3| ) and the relations (|l4|), the exciton Schrodinger equation takes the form 

£n c (fee) + £n h (4) - E? (Q)) ^n h ® + £ / ^ V <^ h ^ ^ ^>'„ = , (15) 

where E^{Q) denotes the energy dispersion of the exciton state. The interaction 

is the in-plane 2D Fourier transform of the 3D Coulomb potential modified due to the confinement in z-dircction. 
The latter is expressed through form factors [with Eq. (14)] 



^L(M') = £ dz e dz h e-\ k - k '\^-^CJze)e: J l {zh)d K {ze)Q J % , (z h ) . (17) 

n e n h mj J J 

The above integrals are calculated analytically since the subband states obtained with the transfer matrix method 
are combinations of exponential and trigonometric functions. 
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The integrable singularity of the Coulomb potential ( |l6| ) at k = k' is taken care of by adding and subtracting in 
Eq. ( |l5| ) the analytically integrable term 

C(k, k') = — ^-[ -r-^ tt-TT^ I • (18) 

1 ^ e 2vr \\k-k'\ rnax(k, k')) ) V ' 

This gives a smooth "corrected" potential (^) of small absolute magnitude. 

To take benefit of the axial approximation, the exciton envelope and the form factors are expanded into 2D angular 
momentum eigenstates exp(i£8). The angular momentum I is chosen for every subband combination such that 1 = 
corresponds to the respective s-like exciton |l3| at Q = 0. 

The resulting set of coupled one-dimensional integral equations is solved numerically for various values of the 
COM momentum Q. Results will be presented in section |y|. 



B. Solution in real space 

We demonstrated in Ref . [jl9| that calculations of the multiband-exciton groundstate dispersion are also feasible with 
a finite-differences scheme in real space. This method is conceptually simple: the Schrodinger equation corresponding 
to Eqs. (0)- leads to a system of four (number of spin components) coupled partial differential equations in 
the four dimensional space (p,z e ,Zh). The resulting eigenvalue problem involves a large sparse complex Hermitian 
matrix with a substantial number, 44, of nonzero off-diagonals. In contrast to the fc-space approach, the method can, 
in principle, be applied effortlessly to very general heterostructures, like quantum wires and quantum dots. The main 
drawback is the need for huge amounts of computer memory. Indeed, the dimension of the matrix to be diagonalized 
scales with the fourth power of the number of grid points per spatial dimension. The most dense grid we used led to 
a matrix of dimension Nn = 4 x 31 x 31 x 61 x 61 = 1.4 • 10 7 (4 is the number of spin components). 

We used the ARPACK |2£| package to calculate a few eigenvalues and eigenvectors at the lower end of the spectrum. 
ARPACK is an efficient implementation |2i| of the Implicitly Restarted Arnoldi Method that can be viewed as a 
synthesis of the Arnoldi/Lanczos process with the Implicitly Shifted QR algorithm. Storage of the (nonzero) matrix 
elements is not required, only a matrix-vector multiplication utility is needed. Multiple eigenvalues, as they occur in 
our case, offer no additional problems. However, if one needs all the members of a multiplet the iteration subspace has 
to be chosen large enough. We find that an iteration subspace of five times the number of the requested eigenvalues 
(rather than the proposed factor of two (2^]) is usually sufficient. This matter was of no concern for our problem, 
since using the symmetry considerations, Eq. (|i"o|), we can derive from a calculated state also the second one of the 
doublet. For the largest matrices, we used a factor of three as a compromise between memory demand and CPU time 
usage. 

The matrix resulting from the discretization is highly structured. For minimizing the memory costs and still making 
full use of the vector registers, we construct the matrix-vector product using auxiliary, much smaller matrices. 

Since memory is critical, it is crucial for any real-space approach to optimize the convergence of the relevant 
quantities with the mesh size. On the one hand, we optimize the COM transformation, as will be discussed in the 
next section, thereby improving the handling of the kinetic terms. On the other hand, we use a groundstate-adapted 
discretization of the Coulomb potential |30| ] which is discussed in more detail in the Appendix. The idea behind 
this approach is to extract the discretized interaction from a reference system that has the same interaction but a 
simple kinetic term, and whose groundstate is known analytically. If the groundstate of the reference system is similar 
enough to the one of the real system, good convergence is expected. 

Calculations on a parallel-vector machine of type CRAYJ932 reached performances of 140MFlops/CPU, the peak 
performance of the limiting BLAS routine being 185MFlops/CPU. 



III. OPTIMIZED CENTER-OF-MASS TRANSFORMATION 

We return to the ambiguity in the COM transformation (^) which is expressed in the freedom to choose /3. The 
relevance of f3 for accelerating numerical convergence in dispersion calculations was realized quite early for bulk 
excitons in where a scalar (3 was optimized in a trial and error procedure. However, there has been no other 
algorithm to take advantage of this freedom until recently ||i"9fl . Before that, there have been just two publications 
where numerical multiband exciton dispersion in quantum wells were calculated: in Ref. jlTj f3 = 1 (in the parabolic 
case m e = oo) was taken in order for the form factors (|T^ ) to be independent of Q, and in Ref. Jlq| no particular 
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choice or handling of (3 is mentioned. In analytic expressions, usually the symmetric (in the parabolic case m e = m/j) 
value (3 = 1/2 is taken pjl. 

The effect of the (3 choice becomes clear when one evaluates Eq. (Q) for two different values (3, f3' = f3 + 8(3 giving 
R' = R + (5(3) p, k' = k — (8(3) Q. Clearly, (3 moves artificially part of the plane wave of the COM motion into the 
relative part of the exciton @ or, equivalently, it shifts the relative part of the wavefunction in /c-space. A good 
choice of f3, as in the parabolic case (||), keeps the relative part of the exciton in real space as smooth as possible or, 
equivalently, pins the relative part of the wavefunction in /c-space to the origin. This situation is illustrated in Fig. |l| 
where we plot the envelope of the HH\C\ — Is exciton in the single subband approximation using the symmetric 
value of (3 = 1/2. This value of (3 is indeed not optimal, as the large shift demonstrates. Diamonds mark where the 
position of the origin would be for other values of (3. 
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FIG. 1. Logarithmic contour plots of the squared HH\C\ — Is exciton envelope in fc-space for a 5 nm GaAs/Alo.3Gao.7As 
QW at Q x — 0.5 nm" 1 plotted in a coordinate system with the not optimized value (3 — 1/2. Diamonds mark, from left to 
right, the origin of the shifted coordinate systems [at (1/2 — (3)Q X ] for (3 = 0, flo, 1/2, 1, with (3a evaluated from Eq. (p4|). The 
Q — envelope is plotted for comparison (dashed). 



We introduced in Ref. |L9| a quasi-analytical method for determining the optimal choice of the scalar (3, which 
we briefly summarize here. It is motivated by the fact that in the parabolic case the correct COM transformation 
decouples the relative motion and COM motion completely. A full decoupling is not possible for non-parabolic 
dispersions. We looked for a choice of (3 that decouples "as much as possible" . To quantify this, we inserted in Eq. (jl]) 
the general /3-dependent COM transformation (Q), separated the Q-dependent terms from the rest 



H 



H(°) +H^((3)Q + H^((3) Q 2 



(19) 



and viewed these as a perturbation of the Q — exciton. Taking into account the inversion symmetry of the Brillouin 
zone, the kinetic mass of the groundstate exciton g is given in second order perturbation theory by 



2Mf 



= (g\H^(f3)\g) + 



E 



(a\HW((3)\g)[ 
25f(0)-££(0) 



(20) 



However, the exciton mass must not depend on (3. Maximizing the first order contribution in (|2(J) and minimizing 
this way the strictly positive contribution of the higher states to the groundstate mass leads to the analytical result 



(3 = (g\H^\g) /(g\HP+HP\g) 



(21) 



-ff^ are simply the material-dependent coefficients of the (3 2 Q 2 terms when inserting ( p^ ) into the A-space represen- 
tation of the kinetic energies in Eq. (||) and (||), respectively fljj| . 

The explicit form of Eq. (p0|) with the contributions from the higher states dropped and (3 from Eq. (|2l]) suggests 



to define COM-related, effective masses 



for electron and hole 



1/mJ 



(g\H;l\g) satisfying M, 



v 



(22) 



G 



Numerical results show that the masses obtained from ( |22| ) tend to be too small. Nevertheless, the obtained values for 
(3 in [ ^9pl}| were quite reasonable because of the much heavier hole mass. If one actually calculates the contributions 



of the higher exciton states to in (EOT) (which are dropped in (22)), one finds that the only important correction 



comes from the coupling to the LH\C\ — ls-like state. Taking in Eq. (|20|) this single correction term into account 
gives practically the exact curvature of the exciton groundstatc dispersion at Q = 0. 

The above procedure is not the best for determining the optimal value of (3 as the importance of the coupling to 
higher states demonstrates. It was inspired by the solution of the exciton problem in real space. Let us now look 
at the form of the exciton wavefunction for Q ^ in the subband expansion (|l3"|). The Q~dependence enters the 
wavefunction: (i) the need to appropriately combine the subband states to get the right Q and (ii) through the need 
for the envelope to adjust for the anharmonicities in the dispersions. In the perturbation approach described above, 
we tried to find a COM transformation that keeps for small Q values the entire wavefunction unchanged as much as 
possible. However, once the one-particle problem is solved, the Q-dependence due to the appropriate combination 
of the subband states (i) is explicitly known. Therefore, a better Ansatz for the wavefunction would be to find a 
COM transformation that keeps the envelopes as much unchanged as possible; that is Eq. (113) with 



( 23 ) 

The minimization of the energy with respect to f3 can be done analytically in the limit Q — > 0. We easily obtain the 
optimized value 

E Jdk\^ nh (k)\ 2 (Q-V) 2 £n h (k) 

A) = • ( 24 ) 

E jdkl^nfnM | 2 (Q • V) 2 (£„ e (fe) +S nh (k)) 
n e n h 

This expression accounts also for the dependence of (3q on the direction Q of the COM momentum in the case 
of warped valence bands. We took into account the inversion symmetry of the Brillouin zone and assumed that the 
Coulomb potential (|l6|), ([Tt]) can be approximated as a function of the momentum transfer only, V(k, k') ~ V(k — k'), 
i.e., we neglected any Q-dependence of the (Coulomb) potential energy of the groundstatc. This is expected to be 
a good approximation since the Coulomb energy depends solely on the charge distribution, which should not be 
affected significantly by the in-plane motion. Indeed, it has been estimated in Ref. [jllj that the error introduced 
by neglecting in Eq. (|l7|) the fc-dependence of the hole envelopes is about 5%. Our assumption should lead to even 
smaller deviations. 

Fig. [j] demonstrates the quality of the expression (|4|). It shows the groundstate envelope at a rather large value 
of Q = 0.5 nm -1 , even though Eq. ( |2~4| ) was obtained in the limit Q — ► 0. It is particularly impressive that the 
respective origin for the optimized choice /3q from Eq. ( p4[ ) lies even a bit to the left of the maximum of the envelope. 
This accounts for the p-component that deforms the originally radially symmetric Q — envelope; using an angular 
momentum decomposition of the envelope a minimum number of components would be needed. The slight deformation 
(these are logarithmic plots) for large values of Q is due to the anharmonicity of the one-particle dispersions. 

The importance of a suitable choice of the COM transformation for the numerical convergence is illustrated in 
Fig. H for the dispersion of the HH\C\ — Is exciton of a 5 nm QW. This has been calculated in fc-space for various 
values of j3 with the same basis (HHxCx, LHiCi, I = 0, ±1, ±2). The further the used (3 lies from the optimal value 
Pa (Po — 0-23 in this case) the worse the results are. We did also calculations where for each Q value an optimal value 
of (3 was obtained by numerical variation. We observed deviations from (3q less than 1% near Q = and not larger 
than 10% at Q = 0.5 nm -1 even for the widest well (20 nm). At large Q the /3 = (in the parabolic case rrih = oo) 
curve gives slightly better results than (3q since the HH\ subband dispersion gets more flat after the avoided crossing 
with the LH\ subband, but it gives considerably worse results at small Q. 
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FIG. 2. Dispersion of the groundstate exciton of a 5 nm wide QW calculated with the same basis but for various values of 
P. The exciton continuum edge (dotted), Eq. is given for comparison. 



IV. A SIMPLE ANALYTICAL FORMULA FOR THE AVERAGE EXCITON KINETIC MASS 

The analytical variation that led to Eq. (|24|) gives the groundstate energy up to terms quadratic in Q. The 
corresponding groundstate kinetic mass Mg has again the form 

Mf =m?+mX (25) 
with the COM-related effective masses for electron and hole defined as 

l Kh = ^Y. fdk\^ nh (k)\ 2 (Q-V) 2 £ e , h (k). (26) 

With these masses, the expression for /3o, Eq. (|24|), has the same form as in the parabolic case (|^). Eq. ( p6| ) gives the 
correct results for the free particle case. 

We claim that this simple result will be of considerable practical importance. Eq. (^6|) is physically appealing: it 
leads to a weighted average of the subband dispersions. Further, it is relatively simple to calculate: It requires only 
the approximate knowledge of the Q — exciton envelope and of the involved subband dispersions. The numerical 
calculation of subband dispersions is nowadays an easy task (provided the k ■ p parameters are known). Moreover, 
especially for narrow QW, the envelope of the HHiCi component of the groundstate exciton is to a very good 
approximation similar in shape to the groundstate of the 2D exciton, 



rf?(p) = J^re-°^, v \° { k) = \l 2 ^(l + (a B kfY m . (27) 
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The LHiCi component is quite small, e.g., 5% for the 20 nm QW, and can be safely neglected in this context. 
Therefore, only a good estimate for the effective Bohr radius as is needed to evaluate (Eq). 



V. RESULTS 



We have calculated exciton dispersions both in real and momentum space for GaAs/Alo.3Gao.7As (001) QW of 
various widths. The coupling of heavy and light holes was fully incorporated. The values of the material parameters 
7ij 72 j 73, m e were taken by linear interpolation from the GaAs and AlAs values; the offset ratio was V v /V c = 0.68/0.32 
and the band gap in meV was taken as E g {x) = 1519 + 1040 x + 470 x 2 , x being the Al content |32| . For the dielectric 
constant, we adopted e = 12 for both well and barrier material. 
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A. Subband expansion 



The nomenclature is as follows: the exciton in the subband expansion has various n^n e subband components with 
the corresponding envelopes Vnln h (F)i Eq. (fl3|). These envelopes have in the axial approximation at Q = a definite 
angular momentum I and will be denoted by Is, 2s, 2p±, 3d± and so on. Each exciton state at finite Q will be named 
according to the main subband component of the corresponding state at Q = 0. That is, speaking of the HH\C\ — Is 
exciton means that at Q = its main subband component is the HH\C\ one with an Is envelope. Similar to the 
single-particle hole subband states, which can change their heavy- or light-hole character away from the T-point, the 
envelope of the main subband component or even the main subband component itself can change with increasing Q. 
To denote the main subband component of a state at a given value of Q, we will speak of the character of the state at 
this Q. For example, the HHiCi — 2p + exciton has a HH\C\ — 2p + character at Q = and a LH\C\ — Is character 

O»o. 

The exciton dispersions in fc-space are calculated as follows: For each QW, we first calculate the exciton spectrum 
at Q = 0. Subsequently, a 2D ls-exciton groundstate function ( p7|) is fitted to the HH\C\ envelope. For the wider 
QW, also a two-dimensional 3d-exciton function is fitted to the LH\C\ envelope. This fit is used, instead of the 
numerical envelope, to evaluate the optimized COM transformation (|2^) because it allows to take advantage of the 
analytically known derivatives of the fit function. The so calculated value of (3q is used for the Q ^ calculations. 

In Fig. y, we display the envelopes of the components of the groundstate exciton in the subband expansion (|l3| ) 
for a 15 nm wide QW at Q = 0. The coupling of the higher subbands is rather small, less than 3% for the LH\C\ 
component and even less for the others. The parity and spin selection rules for the Coulomb coupling of the subbands 
at the T-point in symmetric QW are obeyed, e.g., the admixture of the LH\C\ state vanishes at the T-point since 
the Coulomb potential is spin-diagonal and the HHi and LHi subband states are pure heavy- and light-hole states, 
respectively. The HH\C\ envelope is very well approximated by a 2D ls-exciton function. Deviations are mainly 
located at the vicinity of the HH\ — Lti\ avoided crossing of the hole-subband dispersions (here, k ac ~ 0.13 nm -1 ). 
The total in-plane probability distribution follows the form of Eq. ( |27| ) even better than the HH\C\ envelope alone; 
the coupling to the higher subbands allows the exciton to relax further. This, again, supports the notion that the 
subband mixing has little influence on the charge distribution. 




FIG. 3. On the right: Envelopes of the subband components of the groundstate exciton calculated in fc-space for a 15 nm 
wide QW at Q = 0. On the left: HH\C\ envelope (x), fitted 2D ls-exciton wavefunction (line), and square root of the total 
probability distribution [0(fc)| 2 = J2 n n (p ( ^f„ h (k) 2 (circles). 

The calculations in fc-space presented in the following take into account only the two lowest hole subbands (n^rij, = 
HH\Gx, LH\C\). Inclusion of higher subbands does not enhance the binding energy of the groundstate exciton 
considerably. For the angular decomposition of the envelope components only the s, p±, d+ (£ = 0, ±1, 2) components 
for the HH\C\ and the s, p±, d- (I — 0, ±1, —2) components for the LH\C\ were considered. Due to the optimized 
choice of the COM coordinate system, Eq. (p4|), these few angular momentum components are sufficient to describe 
the HHiC\ — Is and LH\C\ — Is dispersions excellently over the whole range of COM momentum values considered, 
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Q < 0.5 nm 1 : The p components account mainly for the deformation of the envelope, and the d components take 
care of the Coulomb coupling to higher states. 






FIG. 4. Exciton dispersions, groundstate (thick) and some excited bound and continuum states (thin), calculated in fc-space 
for a 15 nm wide QW. The light (dark) gray shaded area mark the HHiCi (LHiC'i) exciton continuum. A parabola for the 
average groundstate exciton mass Mf after Eqs. is plotted (dashed). An arrow marks the position of 1/clb (as the 

Bohr radius). 



In Fig. we show the dispersion calculated in fc-space of the first bound exciton states as well as some of the 
continuum states in a 15 nm wide QW. Zero of energy is the onset of the HH\C\ continuum at Q = 0. The Bohr 
radius as is the one obtained from the same fit of a 2D ls-exciton to the envelope of the HH\C\ component at Q = 
which was used to calculate the optimized /3q. The groundstate is HHiCi — Is. The next state is HH\C\ — 2p + which 
turns into LH\C\ — Is character at the avoided crossing. The LH\C\ — Is exciton is at Q = the fourth excited 
state and shows a substantial mixing with the HHiCi — 3d+ exciton. The HH\C\ — 2p + exciton has at Q = a 
slightly lower energy than the HH\C\ — 2p^ (third state at Q — 0) because it couples to the LH\C\ — Is outside of 
the F-point. 

10 



5 



\ ° 

E 
in 

-5 



-10 



0.0 0.2 0.4 0.0 0.2 0.4 0.0 0.2 0.4 

Q ( nm" 1 ) Q ( nm" 1 ) Q ( nm" 1 ) 

FIG. 5. (a), (b): Comparison of the dispersion of HH\C\ — Is and LH\C\ — Is (20 nm QW only) excitons (thick) with the 
appropriately shifted respective exciton continuum edge (dashed) and the underlying hole-subband dispersion (dot-dashed) for 
two GaAs/Alo.3Gao.7As QW. (c): Comparison of our results (black) with those of Ref. (gray) for the 5 nm QW (shifted 
to match at Q = 0). The exact _H"Hi-subband dispersion for the parameters of |llj] is plotted, too (diamonds). 

As was reported in earlier work [ft7| , [18| , the exciton COM dispersions are highly non-parabolic, much like the hole 
subband dispersions. This is not surprising, since the dispersion of the conduction subband is parabolic and the hole 
mass is much larger than the electron mass. Furthermore, it has been claimed in Ref. that the exciton dispersion 
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follows, in a good approximation, the hole subband dispersion. Although this is certainly true in the present case due 
to the parabolic electron dispersion and the small electron to hole mass ratio, the exciton dispersion is in principle a 
iwo-particle quantity. In fact, the dispersion of the groundstate exciton follows even more closely, in the studied cases 
within 1 meV, the electron-hole-pair continuum edge £ n< , nh (Q)- The latter is defined by 

£n c n h (Q) = min \E n<l (k e ) + £ nh (k h ) \ , (28) 
k e +k h =Q <■ > 

and represents the minimal kinetic energy of a free electron-hole pair for a given subband combination n e rih and 
a given Q. In the independent-subband approximation, this coincides with the exciton continuum edge. Fig. |5|(a,b) 
directly compare for a narrow and a wide QW the exciton groundstate dispersion and the appropriately shifted 
exciton continuum edge. Also shown is the respective hole dispersion. The latter lies always above the shifted exciton 
continuum edge since Eq. (H|) implies £ ne n h (Q) < £n„(0) + £ nh (Q)- 

The exciton groundstate dispersion is found to lie systematically below the shifted exciton continuum edge, i.e., 
the groundstate exciton binding energy becomes larger away from Q = 0. This can be understood based on the fact 
that for Q ^ the groundstate exciton, Eqs. (13 14|), is built from hole subband states around (1 — Po)Q- Due to 



the flatter subband dispersion around this point, the hole mass gets larger and the wavefunction can better adjust 
to the potential. The LH\C\ — Is exciton in the 20 nm wide QW, Fig. ||b, is well separated from the spectrum of 
the HH\C\ exciton. We observe that at small Q the LH\C\ — Is dispersion lies above the respective shifted exciton 
continuum edge. Indeed, the LH\ subband shows a negative mass at the T-point which becomes positive near the 
avoided crossing. Hence, the respective exciton has to pay with extra kinetic energy in order to achieve a small 
COM momentum, and its binding energy is decreased. For larger Q values it gains again some binding energy. 
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FIG. 6. Enhancement of the exciton binding energy Et(Q) with increasing center of mass momentum Q for 
GaAs/Alo.3Gao.7As QW of various widths. The arrows mark the position of the respective HHi-LHi avoided crossing, from 
the right to the left for growing QW width. 



The enhancement of the groundstate binding energy for Q ^ is particularly large when the exciton is built from 
hole subband states around the avoided crossings (1 — (3q)Q = k a . c .- This is demonstrated in Fig. || for QW of 
various widths. Peaks are seen at the respective location of the HH\-LH\ avoided crossing, marked by arrows. The 
enhancement of the groundstate binding energy with Q is less than 15% and is generally larger for wider QW. 
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FIG. 7. Exciton dispersion for 5 nm GaAs/Alo.3Gao.7As QW. The HHiCi (thick) and LH\C\ (dashed) exciton continuum 
edges are also shown. The background shading indicates how strong is the mixture of the LHiC'i exciton states to the 
numerically calculated eigenstates (black corresponds to 100%). For comparison, also a LH\C\ (dot-dashed) exciton continuum 
edge shifted to the energy of the LHiCi — Is state at Q = has been drawn. 

In Fig. 0, we plot the exciton dispersion of a 5 nm QW. The exciton dispersions, like the hole subband dispersions, 
are less anharmonic in the much narrower QW than for the QW of Fig. |^ because of the larger energy separation of 
the hole subbands. The avoided crossing takes place at rather large k (k = 0.31 nm -1 ). The shading indicates the 
percentage of the contribution of the LH\C\ subband states to the norm of the numerical eigenvectors. The very 
light shading of the bound exciton states confirms the small admixture of the LH\C\ states in the groundstate. The 
shading is darker in the vicinity of the avoided crossing. Small-scale intensity variations are numerical artefacts due 
to the finite fc-space mesh. In the exciton continuum only the shading is shown. The LH\C\ — Is resonance is clearly 
seen starting at approximately 18 meV. The resonance is sharper at the beginning and becomes more diffuse at the 
avoided crossing. It lies slightly above the respective exciton continuum edge. 



B. Real-space calculations 



For effectively two-dimensional structures with translational symmetry like the considered symmetric QW in axial 
approxi mati on and, maybe, for some highly idealized quantum wire structures, the real-space method presented in 
section [IB can not compete with the one in fc-space. However, for realistic one-dimensional structures, like V- 



groove and T-shaped quantum wires, this may be the only feasible approach for calculating the exciton groundstate 
dispersion. This is due to the high number of confining dimensions for the exciton (in quantum wires four, two for each 
particle): an expansion in a problem-adapted basis like the product basis of the one-particle eigenstates, for which 
one expects reasonable convergence, leads to four-dimensional integrals for the Coulomb interaction. An expansion 
in a basis where the Coulomb potential is simple will probably show a very slow convergence with basis size. 

The calculations reported here are mainly to be viewed as tests of the applicability of our real-space approach 
and its generalization to finite-elements discretization. They are primarily compared with results obtained with the 
more established fc-space methods. We will therefore discuss the results of the real-space calculations focusing on the 
convergence properties of the method. Further, wavefunction features are better visualized in real space, in particular, 
the electron-hole correlation in the growth direction. 

For the real-space calculations at finite Q, we used the optimized (3q obtained in the respective fc-space calculations. 
We could have used equally well some other procedure to find the effective Bohr radius as, e.g., a variational one, or 
we could also have fitted a 2D ls-exciton function to the in-plane probability distribution of the previously calculated 
exciton groundstate at Q = 0. 

The Coulomb potential was discretized as described in the Appendix. The 3D ls-exciton in the four dimensional 
space (p, z e ,Zh) with m r e e f = 0.0665 mo, tn r ^ = 0.24 mo was used as reference groundstate. The value for the reference 
hole mass was taken from Fig. |ll|, discussed below, as an average value for the range of QW widths considered. This 
gives a reasonable reference Bohr radius of — 12.2 nm; it is nearly the correct value for the in-plane motion or 
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somewhat larger. In the confinement direction the size of the reference wavefunction is larger than the actual one 
(the exciton is quenched in this direction), too. As discussed in the Appendix, a reference Bohr radius as large as 
or somewhat larger than the actual one gives good convergence. We did test calculations with doubled and 
with halved and found a qualitatively similar behavior as in Fig. [jj] in the Appendix. The integration region 
was 120 nm wide in the p directions and 30 nm (60 nm) wide in the z directions for the 5 nm (20 nm) QW with a 
gridpoints distance of about as/6 in the in-plane directions. 












• 









FIG. 8. Middle panel: Exciton dispersion for a 5 nm wide QW. fc-space results for the groundstate and some excited bound 
and continuum states (lines) and the exciton continuum (shaded) are presented, as well as some real-space results (diamonds) . 
Side panels with logarithmic contour plots of the exciton probability distribution for each spin-component (from top to bottom 
raj = +§7+5,— 5,— f) are shown for characteristic exciton states (thick lines) at Q — (left) and Q x = 0.5 nm - (right 
panels). The displayed area is in each direction 10 nm wide for the z e zn and about 10a_g wide for the p x p y plots, respectively. 



The panel in the middle of Fig. g shows dispersions of the lowest exciton states (diamonds) from the real-space 
calculation for a 5 nm wide QW. The exciton dispersions calculated in momentum space are plotted as full curves for 
comparison. The groundstate binding energy is not yet fully converged: all real-space results were shifted approxi- 
mately 1 meV to lower energies to match the groundstate energies at Q — for both methods. Numerical tests show 
that the gridpoint density in the growth direction (« 0.6 points/nm) is more critical than in the in-plane direction 
(« 0.6 points/nm). The exciton continuum edge lies also 0.5 meV too high, and the stronger confinement due to the 
Coulomb interaction lets us expect for the exciton a larger deviation. Nevertheless, the groundstate relative disper- 
sion, £g(Q) — £g(0), is converged and reproduces the fc-space results very well. We note that in the case of parabolic 
one-particle dispersions, this is an exact property of any numerical exciton dispersions. The dispersions of the excited 
states are not reproduced that well. This is mainly due to their larger spatial extension and smaller energy separation 
from each other compared to the groundstate. 

The panels on the left and on the right in Fig. |^ show logarithmic contour plots of the exciton probability distribution 
for some characteristic states. The probability distribution is either integrated over z e , Zh and displayed in the p plane, 
or integrated over p and displayed in the z e Zh plane for each spin component separately. The numerically obtained 
wavefunctions are a linear combination of the two degenerate solutions with opposite parity, Eq. (^|). They are 
disentangled according to parity P, and only P = — 1 states are displayed. 

The left panel displaying the groundstate exciton at Q = illustrates its HH\C\ — Is character: the main spin 
component is mj = +3/2 and has no nodes. The bulk of the exciton is confined in the QW but there is substantial 
penetration into the barrier, being stronger for the lighter electron (note the logarithmic plot). At Q = 0.5 nm -1 
(right panels), the groundstate exciton has still HH\C\ — Is character, as is seen in the fc-space calculations, Fig. 3. 
This does not contradict the strong mixture of heavy- and light-hole bulk states seen in the lower right panel of Fig. |. 
Indeed, the exciton is built from hole subband states near (1 — (3q)Q. This point lies past the HH\ — LH\ avoided 
crossing. Hence, the HH\ subband states near this point are a strong mixture of heavy- and light-hole bulk bandedge 
states. The stronger penetration into the barrier of the light-hole component is again related to its smaller mass. 

An interesting feature is the larger confinement of the mj = +| component at Q — 0.5 nm -1 compared to Q = 0. 
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This is a consequence of the enhanced exciton binding energy. Altogether, the plots demonstrate that the total 
charge distribution in not altered much with increasing COM momentum; the somewhat stronger confinement of the 
heavy hole is at least partly canceled by the larger penetration into the barrier of the light hole. Although here not 
clearly resolved, the in-plane plots at Q = 0.5 nm -1 show the slight deformation of the originally radially symmetric 
wavefunction that was seen in Fig. [lj. Again these deformations partly cancel each other in the sum over the spin 
components, and the in-plane charge distribution remains mainly symmetric. This is more clearly seen for the first 
excited state at Q = 0.5 nm -1 which has H H\Ci-2p y character (upper right panel in Fig. ||). It is the HH±Ci — 2p+ 
exciton, which is, again, slightly lower in energy than the HH\C\ — 2p_ exciton at Q = 0. At large Q the character 
of the HH±Cx envelope changes from 2p± to 2p VtX . 




z 

e 



FIG. 9. Groundstate exciton probability density in the z e Zh-plane for the main spin-component (raj = +|)atQ = for a 
5 nm wide QW, (a) integrated over the in-plane relative coordinate p (same as lower left corner of Fig. ^) and (b) at p — 0. 
The background lines mark the position of the QW interfaces. 

One does not expect strong electron-hole correlation in the growth direction for the 5 nm QW, which is considerably 
narrower than the exciton Bohr radius. Here, the confining potentials arc on average much stronger than the Coulomb 
potential and the wavefunction can not relax in this direction. Indeed, in Fig. ^(a) the probability density integrated 
over p does not show much correlation: the contour lines are not elongated along the diagonal, z e = Zh- However, some 
correlation exists as the plot of the cut at p = in Fig. ^(b) demonstrates. This weak correlation is not included in our 
fc-space calculations where only one conduction band was used. Its envelope does not depend on k and consequently 
the z e coordinate can be separated, Eq. (|l3|). 




FIG. 10. Logarithmic contour plots of the exciton probability distribution in the 2 e Zh-plane for each spin-component (in 
each panel from top to bottom mj — +|, + |, — |, — |) for the HH\C\ — Is and the LH\C\ — Is excitons at some characteristic 
Q-values for a 20 nm QW. The displayed region is in each direction twice the QW width. The triangles (diamonds) show results 
of the real-space calculations for a dense (less dense) mesh. Lines and shading of the background panel are the fc-space results 
of Fig. El We do not show the HH2C1 continuum that lies partly in the displayed region. 
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In Fig. [lO[ results are displayed for a 20 nm wide QW. Two set of points are shown for two different mesh sizes. 
The energies of the groundstate exciton are almost converged for the more dense mesh; the deviation from the results 
of the fc-space calculations is only 0.2 meV at Q — 0. 

The small, but noticeable mj = +1/2 component of the groundstate exciton HH\C\ — Is exciton in the lower 
left panel corresponds to the substantial admixing of the LH\C\ — d- exciton seen already for the 15 nm QW in 
Fig. y. The admixing is larger for wider QW due to the smaller energy separation of the respective subbands. At 
Q = 0.5 nm -1 , it has mainly bulk light-hole character (mj = —1/2). This, again, does not contradict the HH\C\ — Is 
character since the HHi subband has approximately 60% light-hole character beyond the avoided crossing with the 
LHi subband. The first excited state at Q = 0.5 nm -1 (LH±Ci — Is exciton) has LH\C\ — Is character even though an 
additional node is seen in the z e Zh probability distribution of the main spin component (mj = +1/2). The envelopes 
of the single-particle subband states at large enough in-plane momentum show more nodes than at the T-point due 
to the coupling of in-plane and growth directions in the Luttinger Hamiltonian. This is one of the reasons why the 
expansion Eqs. (|l3|,[l4|) gives very good results with just two subbands, while an expansions in the subband states at 
the r-point |18| needs more subbands for the same accuracy. 

The probability distribution plots illustrate the almost vanishing penetration into the barriers for the wide QW, in 
contrast to the narrower QW of Fig. |[ All panels show a clear orientation of the contour lines towards the diagonal, 
z e = Zh. This demonstrates the considerable electron- hole correlation in the growth direction for QW wider than one 
Bohr radius. Still, the stronger correlation has little impact on the energies. Recall that in perturbation theory the 
first order correction to the wavefunction gives only a second order correction to the energy. This justifies the usual 
factoring out of the dependence on the growth direction for the much lighter electron. 

This, however, is not the case for the much heavier hole. It is a well known fact, that the Coulomb coupling of the 
hole subbands is considerable. Indeed, for the 20 nm wide QW, neglecting the Coulomb coupling of the HHiC\- to 
the L-ffiCi-excitons leads to an error in the groundstate binding energy larger than 10%. That is, the correlation of 
in-plane and confinement directions for the hole is substantial. 

For the 20 nm QW, a single calculation at Q = 0.5 nm -1 was done with /3 = 1 in order to check the relevance 
of this parameter for the numerical accuracy in the real-space calculations, too. Indeed, the respective groundstate 
energy lies very far ( 43 meV ! ) above the correct value. 



C. Average exciton groundstate mass 



In the process of determining the optimal choice of the COM coordinate system, we derived in section 111 an 
expression for the kinetic mass of the groundstate exciton, Eqs. ( p5| , p6| ). Two assumptions were essential: First, the 
Coulomb potential, i.e., the form factors, are a function of the in-plane momentum transfer, and, second, the Ansatz 
(p3|) is valid. The numerical results of the previous sections support these assumptions. 

Together with the exciton dispersions in Fig. ^, we displayed for the groundstate exciton a parabola with the exciton 



mass Mf of section pi This mass is obviously not the one determined by the curvature of the exciton groundstate 



dispersions at Q = 0. It is rather an average of the curvature of the grou ndst ate dispersion in a region of size l/a# 



Indeed, this is implied by Eqs. (|25|,[26|) and the observation made in section V A that the groundstate dispersion follows 
closely the respective exciton continuum edge. 

This exciton COM momentum region is the one important for exciton localization due to well width fluctuations, 
interface roughness or alloy fluctuations in QW Indeed, the exciton averages over smaller scale fluctuations due 
to its finite size and feels an effective disorder potential that is spatially correlated over the Bohr radius a^. 
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FIG. 11. Exciton average kinetic mass in various appro ximations and Bohr radius as versus well width L for 
GaAs/Alo.3Gao.7As QW. The figure is discussed in section |V Q 

The dependence of various expressions for the groundstate exciton mass on the QW width L is displayed in Fig. [Tl| . 
The conduction band mass was taken material-independent m f ; w = 0.0665, the shown L-dependence comes solely 
from the valence band. Displayed are the masses obtained by: (i) describing the hole in the diagonal Luttinger 
approximation 1/m/i = P w (7™ + 7^) + Pb (7J + 72), where P W) b denote the probability that the hole is in the well 
and barrier material respectively (long dashed) , (ii) taking as hole mass the subband curvature at the T-point that is 
known analytically [^3|l (dashed), (iii) using our semi-analytical expression ( p5|j26|) for the mass with a 2D ls-exciton 
function fitted to the envelope of the HH\C\ component (diamonds), and (iv) the average of the curvature of the 
numerical exciton dispersion weighted with the same function as (iii) (circles, the line is a guide to the eye). 

The top curve displays the values for the Bohr radius that we have used for the calculation of (iii) and (iv). Arrows 
at either sides of the lower part of Fig. [ll] mark the mass in the diagonal approximation in the well and barrier bulk 
materials. 

Fig. [ll] demonstrates the failure of the diagonal Luttinger approximation (i) to describe even the HH\ subband 
curvatures at the T-point due to the degeneracy of heavy- and light-hole bands in the unstrained bulk. However, 
even the correct single-particle subband curvatures at the T-point (ii) fail to describe accurately the curvature of the 
groundstate exciton dispersion at Q = 0. This is mainly due to the finite extension of the exciton in fc-space that 
implies an averaging of the subband dispersions over a region of approximately 1 /as near the T-point and partly to 
the Coulomb coupling to higher subbands. Both effects tend to make the groundstate exciton heavier. The mass 
derived from the curvature of the groundstate exciton at Q — (not shown) lies between curves (ii) and (iii) . 

The numerically obtained "best" mass values (iv) show a quantitatively and qualitatively different behavior from 
curves (i) and (ii). For very narrow QW the subbands become flatter at the T-point because of the larger penetration 
into the barriers where the exciton becomes again heavier, as in models (i), (ii). But, for large L the region of the 
H H1-LH1 avoided crossing comes to a distance of approximately I /as to the T-point and the exciton, averaging over 
the flatter subband dispersion, becomes heavier. 

The quality of our semi-analytical result for the average exciton mass (iii) has to be judged according to its deviation 
from the numerical average (iv). The non-monotonous behavior is clearly seen for the mass (iii) obtained using only 
the HHi subband dispersion and the fitted Bohr radius. The mass values (iii) are somewhat smaller (<10%) than the 
ones of curve (iv) . This is due to the enhancement of the binding energy for Q 7^ that yields larger average masses 
(iv) than one would expect based on the one-particle subband dispersions and the Q — groundstate exciton. Our 
semi-analytical average groundstate exciton mass (iii) reaches a minimum approximately at the QW width where the 
maximum binding energy is reached. 

The small differences between curves (iii) and (iv) demonstrate the quality of our expression (HJH) for the average 
exciton groundstate mass. For this reasonable and easy-to-use mass expression, only a good estimate for the in-plane 
Bohr radius and the dispersion of the involved single-particle subbands is needed. 

In the paper by Triques and Brum fjs|] , average exciton effective masses were calculated that are relevant to the 
formation process of excitons in two different scenarios. These masses are defined by parabolic fits within a relevant 
energy range of: (a) 5 meV, approximately half the exciton binding energy, in the case that the particles first relax 
and then form an exciton of kinetic energy lower than the binding energy and (b) 36 meV in the case that the exciton 
is formed very fast and relaxes initially via optical phonon emission, reaching energies below 36 meV, the energy 
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of the GaAs-LO phonon. These energy ranges translate to Q values in general much larger than 1/ob- Therefore, 
these average masses should be larger than those of the present work. However, the values published in |1| for the 
5 meV mass (0.2mo < < 0.3mo for L < 10 nm) are smaller than ours for narrow QW. Our average masses of 
scenario (a) are not falling below 0.3mo and show a smooth minimum for a QW of width somewhere between 2 nm 
and 5 nm. This difference can be traced back mainly to the inefficiency for narrow QW of their method involving an 
expansion in the subband states at the T-point, as already remarked by the authors themselves. Indeed, for narrow 
QW the few states at the T-point can not provide the needed flexibility to simulate states far away from the T-point. 
This is substantiated by the lower left panel in Fig. |§| showing the groundstate exciton at Q = 0.5 nm _1 for the 5 nm 
QW. The z e Zh-plots show a substantial mj = +^ spin component. In an expansion of the Zh dependence in the hole 
subband states at the T-point this component would need a LH 2 subband to be described efficiently. However, for 
the 5 nm QW only the HH\, LH\, HH2 subbands exist below the top of the barrier. 

In Fig. ||(c), we compare the results of Ref. jll| with ours for the 5 nm QW. For this narrow QW, their exciton 
and HH\ subband dispersions are much steeper than ours. The discrepancy is not due to the different parameters as 
a comparison with the subband dispersion calculated exactly with the transfer matrix method for the parameters of 
Ref. [[18| demonstrates (diamonds). 

We remarked already, that different parts of the non-parabolic groundstate exciton dispersion are relevant for 
different physical process, and we gave exciton localization [QJJ] and cooling of a non-thermal exciton population |l§|l 
as examples. We would like to discuss two more possible experimental consequences of anharmonicities using FigTH 
as illustration. First, an exciton population loosing energy by acoustic phonon emission at low temperatures could 
experience a bottle-neck effect, i.e., an increased population of the fc-space region near 0.15 nm -1 . Scattering closer 
towards the T-point will be suppressed by the decreased density of final states. A second, more directly observable 
consequence is the temperature dependence of the exciton lifetime. Evaluation of the latter along the lines of Ref. |34| ] 
with the groundstate dispersion of Fig. ^ yields a super-linear increase of the exciton lifetime with temperature: The 
slope increases in the range from 5 to 20 K by a factor 2.2 (not shown). 

VI. CONCLUDING REMARKS 

In summary, we performed k ■ p multiband exciton dispersion calculations of high accuracy even for very large 
COM momentum and narrow QW using the well known expansion in subband states. For the high quality of the 
numerical results, an optimized COM transformation was essential. The optimized COM transformation is related 
to an average exciton kinetic mass for which a simple semi-analytical expression is derived. This eliminates for many 
practical purposes the need to actually calculate the groundstate dispersion. The groundstate exciton dispersion is 
found to follow closely the respective exciton continuum edge. 

In addition, we demonstrated that multiband groundstate exciton dispersion calculations arc feasible in real space 
using a finite-differences scheme. Besides the essential optimized COM transformation, a convenient groundstate- 
adapted discretization of the Coulomb potential enhances numerical accuracy. This method promises to give results 
for systems where the first one is not practicable. Its generalization for finite-element schemes is straightforward. 
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APPENDIX: DISCRETIZATION OF THE COULOMB POTENTIAL 

In order to attain a simple discretization for the interaction of some Hamilton operator, Glutsch, Chemla, and 
Bechstedt [[30| proposed to discretize on the same mesh another operator whose groundstate is analytically known 
with the same interaction but with a simple kinetic term (mass m re ^). We will call this the reference system. For 
illustration, consider a simple one-dimensional system with known groundstate wavefunction g(x) of energy E r g e f . 
The Schrodinger equation of the reference system discretized on a mesh xi — i ■ reads 



h 2 



2m re f Al 



(g(x i+1 ) - 2g(x i )+g(x i -i)) + V(x i )g(x i ) = E™* g( Xi ) (29) 
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and yields the groundstate adapted discretization 



V(x t ) = + 



2m re f Al 



](x i+ i) - 2g(xj) + g(xj-i) 



gfa) + o 



(30) 



This procedure is very simple, easy to implement, cheap to calculate, and gives for the reference system always the 
correct groundstate energy, regardless how inappropriate the mesh is. In addition, no special handling for potentials 
with integrable singularities is needed. 




1 2 3 4 5 6 7 
Gridpoints within 1 a B0 



1 2 3 

Gridpoints within 1 a M 



FIG. 12. The lowest eigenvalues calculated in real space of the ideal 3D exciton (left) and 2D exciton (right) of Bohr radius 
dsn plotted versus grid density. Discretization of the Coulomb potential as described in the Appendix with reference Bohr 
radius a^: aso (thick), 2 aso (dashed), a.en/2 (long dashed). For the 2D case, also results are shown for the groundstate with 
the potential integrated analytically in every mesh cell, Eq. (^), (thin) and for the groundstate with the "naive discretization" 
(dot-dashed). The former curve is hardly to distinguish from the one for = 2o_bo- 



Using this discretization for a reference state similar enough to the one sought, one can expect good convergence 
with mesh size. In order to check how the dependence of the potential discretization on the reference state influences 
the results, we performed calculations of the ideal 2D and 3D excitons using discretizations of the Coulomb potential 
based on reference excitonic groundstates with various Bohr radii. Fig. [L2| shows the lowest numerical eigenvalues 
as a function of gridpoint density. One can see that the correct estimation of the "unknown" groundstate is not so 
critical: a reference Bohr radius cib within a factor of two from the actual one, asrjj still gives good results for the 
groundstate for reasonable mesh densities. Apparently, especially for the 2D case, it is better to choose as rather 
larger than smaller in order to get good results also for the excited states. 

We also show for the 2D case the results obtained with the potential integrated analytically on each Cartesian mesh 
element, using 



dx dy 



y In [x + \/ x 1 + y 2 ^ + x In ( y + \J x 2 + y 2 ^ 



(31) 



This gives also a very good convergence with the mesh size. The 2D result obtained for the groundstate with the 
potential integrated analytically only at the origin, and taking everywhere else its value at each gridpoint ('naive' 
discretization) is included in Fig. [12L Since both discretizations are the same at the origin the difference does not 
originate from the divergence at this point. The superiority of the discretizations (|30| ) and ( |3l| ) is obvious. We 
therefore expect Eq. (|3(]) to yield good results even with not very dense meshes. 

For the 3D case a similar expression to Eq. (|3l|) for the Coulomb potential integrated analytically on a rectangular 
box can be derived. However, we could not use this result to attain an alternative discretization for the real-space 
QW calculations because the natural mesh (p Xl p y , z e ,Zh) is not Cartesian in the relative coordinate z e — Zh- 



18 



[1] E. Runge and R. Zimmermann, Festkorperprobleme/ Advances in Solid State Physics 38, 251 (1998); R. Zimmermann and 

E. Runge, phys. stat. sol. (a) 164, 511 (1997) 
[2] R. Kumar, A. S. Vengurlekar, S. S. Prabhu, J. Shah, and L. N. Pfeiffer, Phys. Rev. B 54, 4891 (1996); B. Deveaud, F. 

Clerot, N. Roy, K. Satzke, B. Sermage, and D. S. Katzer, Phys. Rev. Lett. 67, 2355 (1991); T. C. Damen, J. Shah, D. Y. 

Oberli, D. S. Chemla, J. E. Cunningham, and J. M. Kuo, Phys. Rev. B 42, 7434 (1990); Ph. Roussignol, C. Delalande, A. 

Vinattieri, L. Carraresi, and M. Colocci, Phys. Rev. B 45, 6965 (1992) 
[3] M. Z. Maialle, E. A. Andrada e Silva, and L. J. Sham, Phys. Rev. B 47, 15776 (1994); A. Vinattieri, J. Shah, T. C. Damen, 

D. S. Kim, L. N. Pfeiffer, M. Z. Maialle and L. J. Sham, Phys. Rev. B 50, 10868 (1994); J. Puis, F. Henneberger, M. Rabe, 

and A. Siarkos, J. of Cryst. Growth 184/185, 787 (1998) 

H. F. Hess, E. Betzig, T. D. Harris, L. N. Pfeiffer, and K. W. West, Science 264, 1740 (1994) 

G. Dresselhaus, J. Phys. Chem. Solids 1, 14 (1956) 
M. Altarelli and N. O. Lipari, Phys. Rev. B 15, 4898 (1977) 

H. Chu, G. D. Sanders, and Y.-C. Chang, Phys. Rev. B 36, 7955 (1987) 

P. Michler, A. Hangleiter, A. Moritz, G. Fuchs, V. Harle, and F. Scholz, Phys. Rev. B 48, 11991 (1993) 
D. A. Broido and L. J. Sham, Phys. Rev. B 31, 888 (1985) 

D. A. Broido and L. J. Sham, Phys. Rev. B 34, 3917 (1986) 
G. D Sanders and Y.-C. Chang, Phys. Rev. B 35, 1300 (1987) 
L. C. Andreani and A. Pasquarello, Phys. Rev. B 42, 8928 (1990) 

C. Y.-P. Chao and S. L. Chuang, Phys. Rev. B 48, 8210 (1992) 
B. Zhu and Y.-C. Chang, Phys. Rev. B 50, 11932 (1994) 
J. Chen, R. Chen, and K. K. Bajaj, Phys. Rev. B 50, 10947 (1994) 
R. Winkler, Phys. Rev. B 51, 14395 (1995) 

S. Jorda, U. Rossler, and D. Broido, Phys. Rev. B 48, 1669 (1993) 
A. L. C. Triques and J. A. Brum, Phys. Rev. B 56, 2094 (1997) 

A. Siarkos, E. Runge, and R. Zimmermann, Annalen der Physik (Leipzig) 7, 523 (1998) 

E. Kapon, D. M. Hwang, and R. Bhat, Phys. Rev. Lett. 63, 430 (1989); O. Stier and D. Bimberg, Phys. Rev. B 55, 7726 
(1997) 

J. M. Luttinger, Phys. Rev. B 102, 1030 (1956) 

B. Zhu and K. Huang, Phys. Rev. B 36, 8102 (1987) 

D. B. Tran Thoai, R. Zimmermann, M. Grundmann, and D. Bimberg, Phys. Rev. B 42, 5906 (1990) 

E. L. Ivchenko, A, Y. Kaminski, and U. Rossler, Phys. Rev. B 54, 5852 (1996); T. Hassenkam, S. Pedersen, K. Baklanov, 

A. Kristensen, C. B. Sorensen, P. E. Lindelof, F. G. Pikus, and G. E. Pikus, Phys. Rev. B 55, 9298 (1996); R. Winkler 
and U. Rossler, Phys. Rev. B 48, 8919 (1993) 

L. C. Andreani, A. Pasquarello, and F. Bassani, Phys. Rev. B 36, 5887 (1987) 
M. A. Kanehisa, Physica 117B & 118B, 275 (1983) 

L. R. Ram-Mohan, K. H. Yoo, and R. L. Aggarwal, Phys. Rev. B 38, 6151 (1988) 

ARPACK and the parallel version PARPACK can be downloaded free of charge at http://www.caam.rice.edu/software 
U. Eisner, V. Mehrmann, F. Milde, R.A. Romer, and M. Schreiber, SIAM J. Sci. Comp 20, 2089 (1999) 
S. Glutsch, D. S. Chemla, and F. Bechstedt, Phys. Rev. B 54, 11592 (1996) 

Unfortunately, in Fig. 1 of Ref. jl9| a different convention was used, with f3 replaced by /3 — 1/2. 

Landolt-Bornstein, New Series, Part IH/2a, edited by O. Madelung, Springer, Berlin, 1982; A. Onton, Festkorperprobleme 
XIII, 59 (1973); E. Hess, I. Topol, K.-R. Schulze, H. Neumann, and K. Unger, phys. stat. sol. (b) 55, 187 (1973) 

B. A. Foreman, Phys. Rev. B 49, 1757 (1994) 

L. C. Andreani, in Confined Electrons and Photons, edited by E. Burstein and C. Weisbuch, Nato ASI Series Vol. 349, 
Plenum, New York, 1995, p. 57-112 



19 



